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Abstract. For more than ten years, solar-like oscillations have been detected and 
frequencies measured for a growing number of stars with various characteristics (e.g. 
different evolutionary stages, effective temperatures, gravities, metal abundances ...). 

Excitation of such oscillations is attributed to turbulent convection and takes place 
in the uppermost part of the convective envelope. Since the pioneering work of Gol- 
dreich & Keely (1977), more sophisticated theoretical models of stochastic excitation 
were developed, which differ from each other both by the way turbulent convection 
is modelled and by the assumed sources of excitation. We review here these different 
models and their underlying approximations and assumptions. 

We emphasize how the computed mode excitation rates crucially depend on the way 
turbulent convection is described but also on the stratification and the metal abundance 
of the upper layers of the star. In turn we will show how the seismic measurements 
collected so far allow us to infer properties of turbulent convection in stars. 



1 Introduction 

Solar p-modes are known to have finite lifetimes (a few days) and very low 
amplitudes (a few cm/s in velocity and a few ppm in intensity). Their finite 
lifetimes result from several complex damping processes that are so far not clearly 
understood. Their excitation is attributed to turbulent convection and takes 
place in the upper-most part of the Sun, whic h is the place o f vigorous and 
turbulent motions. Since the pioneering work of lLightriilll |l952f) . we know that 



a turbulent medium generates incoherent acoustic pressure fluctuations (also 
called acoustic "noise"). A very small fraction of the associated kinetic energy 
goes into to the normal modes of the solar cavity. This small amount of energy 
then is responsible for the small observed amplitudes of the solar acoustic modes 
(p modes). 

In the last decade, solar-like oscillations have been detected in numerous 
stars, in different evolutionary stag es and with different metallicity (see recent re- 
view bv lBedding fc Kieldsenll2007l ). As in the Sun, these oscillations have rather 
small amplitudes and have finite lifetimes. The excitation of such solar-like oscil- 
lations is attributed to turbulent convection and takes place in the outer layers 
of stars having a convective envelope. 

Measuring mode amplitudes and the mode lifetimes permits us to infer V, 
the energy supplied per unit time into the acoustic modes. Deriving V puts con- 
straints on the theoretical models of mode e xcitation by turbule nt convection 
(lLibbrechdll988l) . Towever, as pointed-out bv lBaudin et al. I (|2005l) . even for the 
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Sun, inferring V from the seismic d ata is not a trivial task 
mic data, this is even more difficult (| Samadi et al. 2008) 



For stellar seis- 
We discuss here the 
problems we face in deriving reliable seismic constraints on V '. 

A first attempt to explain th e observed solar five minute oscillations was car- 
ried out by lUnno fc Katol ( 19621 ). They have considered monopole H and d ipoljl 
source terms t hat arise from a n isothermal stratified atmosphere. ISteinl (| 19671 ) 
has generalised iLighthilll (|1952l )'s approach to a stratified atmosphere. He found 
that monopole source terms have a negligible contribution to the noise genera- 
tion compared to the quadrupole source ternH. Among the quadrupole source 
terms, the Reynolds stress was expected to be the major source of acoustic wave 
generation. It was only at the beginning of the 1970's that solar five minutes 
oscillations have been c l early identified as global resonant modes (|Ulrichl 1970 : 
Leibacher fc Steinlll97lt rDeubnerlll975h . A few years later, iGoldreich fc Keelevl 



(jl977bl GK hereafter) have proposed the first theoretical model of stochastic 
excitation of acoustic modes by the Reynolds stres s. Since this pioneering work , 
different improved models have been developed (iDolginoy fc Muslimiov _ 1984 • 
Balmforth|[l992at IGoldreich et al.lll994t ISamadi fc Goupi]||200ll; IChaplin etal 



2005HSamadi et al.ll2003aHBelkacem et al.ll2006bl . 120081 ). These approaches dif- 



fer from each other either in the way turbulent convection is described or by the 
excitation process. 

In the present paper, we briefly review the different main formulations and 
discuss the main as sumptions and appr oximations on which these models are 
based. As shown bv lSamadi et al.l ( 2003a ). the energy supplied per time unit to 
the modes by turbulent convection crucially depends on the way eddies are tem- 
porally correlated. A realistic modeling of the eddy time-correlation at various 
scale lengths then is an important issue, which is discussed in detail here. We 
will also highlight how the mean structure and the chemical composition of the 
upper convective envelope influence the mode driving. Finally, we will summa- 
rize how the seismic measurements obtained so far from the ground allow us to 
distinguish between different dynamical descriptions of turbulent convection. 



2 Mode energy 

We will show below how the energy of a solar-like oscillation is related to the 
driving and damping process. The mode total energy (potential plus kinetic) is 
by definition the quantity: 

E osc (t) = / d 3 xp v 2 osc (r,t) (1) 

where v osc is the mode velocity at the position x, and po the mean density. 

1 A monopole term is associated with a fluctuation of density 

2 A dipole term is associated with a fluctuation of a force 

3 A quadrupole term is associated with a shear 
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Mode damping occurs over a time-scale much longer than that associated 
with the driving. Accordingly, damping and driving can be completely decoupled 
in time. Furthermore, wc assume a constant and linear damping such that 

dv osc (t) 

= -J7«OBc(t) ( 2 ) 



dt 

where 77 is the (constant) damping rate. The time derivative in Eq. ([2]) is per- 
formed over a time scale much larger than the characteristic time over which the 
driving occurs. 

Let V be the amount of energy injected per unit time into a mode by an 
arbitrary source of driving (which acts over a time scale much shorter than 
I/77). According to Eqs. ([1]) and (2}, the variation of E osc with time is given by: 

^(t)=V-2 V E osc (t). (3) 

Solar-like oscillations are known to be stable modes. As a consequence, their 
energy cannot growth on a time scale much longer than the time scales associated 
with the damping and driving process. Accordingly, averaging Eq. (|3|) over a long 
time scale gives: 



*t=w-o, (4) 

where () refer to a time average. From Eqs. ([3]) and we immediately derive: 

E OBC = ^ (5) 

We then clearly see with Eq. ([5]) that a stable mode has its energy (and 
thus its amplitude) controlled by the balance between the driving (V) and the 
damping (77). Then, the major difficulties are to model the processes that are 
at the origin of the driving and the damping. For ease of notation, we will drop 
from now on the symbol () from E osc and V . 



3 Seismic constraints 

As we shall see later, the mode displacement, Sr osc , can be written in terms of 
the adiabatic eigen-displacement £, and an instantaneous amplitude A(t): 

5r osc =^ (A(t)€(r)e- to -* + cc) (6) 

where cc means complex conjugate, w osc is the mode eigenfrequency, and A(t) is 
the instantaneous amplitude resulting from both the driving and the damping. 
Note that, since the normalisation of £ is arbitrary, the actual intrinsic mode 
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amplitude is fixed by the term A(t), which remains to be determined. The mode 
velocity, v osc , is then given by: 

(r, t) = = ~{-ioj osc A(t) £(r) e—* + cc) (7) 

where cc means complex conjugate. Note that we have neglected in Eq. 
the time derivative of A. This is justified since the mode period (27r/oj osc ) is in 
general much shorter than the mode lifetime (~ 1 / 77) 

From Eqs. and (D, we derive the expression for the mean mode energy: 

Ease = J d 3 xp v^' c = i| A \ 2 Iuj osc 2 , (8) 



where 

r-M 



1= / d?xpo£*.£ (9) 







is the mode inertia. For the sake of simplicity, we will from now on only consider 
radial modes. According to Eq. I|7|. the mean-square surface velocity associated 
with a radial mode measured at the radius r/j, is then given by the relation 



v 2 s (r h ) = -\A\ 2 co 2 osc \ Ur h )\ 2 (10) 

where £ r is the radial component of the mode cigenfunction. It is convenient and 
common to define the mode mass as the quantity: 

where 77, is the radius in the atmosphere where the mode is measured in velocity. 
According to Eqs. (8J, (fl7J| . and (JTTJ) , we derive the following relation: 

E osc = M v 2 s (12) 

It should be noticed, that although M. and v s depend on the choice for the 
radius r/i, E osc is by definition intrinsic to the mode (see Eq. |Tj)) and hence is 
independent of r^. 

Using Eqs. (0, and (fT2)) . we finally derive: 

v 2 s (r h ,uj oac ) = -^J^f (13) 

where F — tj/tt is the mode linewidth, and r\ the mode damping rate. From 
Eq. (H2J), one again sees that the mode surface velocity is the result of the 
balance between excitation (V) and the damping (77 = _T7r). However, it also 
depends on the mode mass (M): For a given driving (V) and damping (r), the 
larger the mode mass (or the mode inertia), the smaller the mode velocity. 

When the frequency resolution and the signal-to-noise are high enough, it is 
possible to resolve the mode profile and then to measure both r and the mode 
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height H in the power spectral den sity (generally give n in m 2 /Hz). In that case 
v s is given by the relation (see e.g. lBaudin et al.ll2005T ): 



v 2 s (r h , w osc ) = 7r C obs H r 



(14) 



where the consta nt Cobs takes the ob servational technique and geometrical effects 
into account (see lBaudin et al.ll2005l ). From Eq. (fTB")) and (|14p . one can then infer 
from the observations the mode excitation rates (V) as: 



V(uS) = 2tt M r v 2 = 2ir 2 M C ohs H r 2 



(15) 



Provided that we can measure r and H, it is then possible to constraint V . 
However, we point out that the derivation of V from the observations is also based 
on models since M. is req uired. Furthermore, there is a strong anti-c orrelation 
between H and r (see e.g. IChaplin et al.lll998t IChaplin fc Basull2008T ) that can 
introduce important bias. This anti-correlation vanishes when considering the 
squared mode amplitude, v 2 , since v 2 s oc H T (see Eq. (TT41 V However, V still 
depends on r, which is strongly anti-correlated with H . 

As an alte r nativ e to comparing theoretical results and observational data, 
Chaplin et all lj2005h proposed to derive H from the theoretical excitation rates, 
V, and the observed mode line width, _T, according to the relation: 



H 



V 



(16) 



However, as pointed-out bv lBelkacem et al" ( 2006b ). H strongly depends on the 
observation technique. The quantity C bs H, is less dependent on the observa- 
tional data but still depends on the instrument since different instruments probe 
different layers of the atmosphere (see below). Therefore, one has difficulty to 
compare values of H C Q bs coming from different instruments. 



3.1 Solar seismic constraints 



Baudin et all (|2005l ) have inferred the solar p-mode excitation rates from differ- 
ent instruments, namely GOLF on-board SOH O, the BiSON and GONG net- 
works. As pointed out by iBaudin et al. (|2005l ). the layer (rh) where the mode 
mass is evaluated must be properly estimated to derive correct values of the exci- 
tation rates from Eq. (| 1 5|) . Indeed solar seismic observations in Doppler velocity 
are usually measured from a given spectral line. The layer where the oscillations 
are measured then depends on the height (rh) in the atmosphere where the line 
is formed. Different instruments use different solar lines and then probe differ- 
ent regions of the atmosphere. For instance, the BiSON instruments use the KI 
line whose height of formation is estimated at the optical depth r sa 0.013. The 
optical dept h associated wi th the different spectral lines used in helioseismology 



are given m 



Houdekl (l2006h with associated references. 



Solar p-mode excitation rates, V, derived by lBaudin et al. ( 2005h are shown 
in Fig. Q] (left panel). For v < 3.2 mHz, ;p GONG and p BiSON are consistent with 
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each other, whereas 'p GOLF is systematically smaller than 'p GONG and 7? BlSON ) 
although the discrepancy remains within 1-er. At high frequency, differences be- 
tween the different data sets are more important. This can be partially attributed 
to the choice of the layers where M are evaluated. Indeed, the sensitivity of 
M. to rh is the larger at high frequency. On the other hand, low-frequency mode 
masses are much less sensitive to the choice of rv Accordingly, the discrepancy 
seen at low frequency between GOLF and the other data sets sug gests that the 
absol ute calibration of the GOLF data may not be correct (see iBaudin et all 
2005h . In Fig. CQ we then present V derived from GOLF data after multiplying 
them by a factor in order that they match at low frequency , p GONG and ■p BlSON . 
We find a rather good agreement between •p GOLF and 'p BlSON whereas, at high 
frequency, , p GONG are systematically lower than "p GOLF ore •p BlSON . The resid- 
ual high-frequency discrepancy is likely due to an incorrect determination of 
the layer where the different seismic measurements originate (see a detailled 



discussion in lBaudin et al.ll2005f ). 



3.2 Stellar seismic constraints 

Seismic observations in Doppler velocity of solar-like pulsators are performed 
using spectrographs dedicated to stellar seismic measurements (e.g. UCLES, 
UVES, HARPS). Such spectrographs use a large number of spectral lines in 
order to reach a high enough signal-to-noise ratio. In the case of stellar seismic 
measurements, it is then more difficult than for heliose ismic observat i ons t o 
estimate the effective height r^. As discussed in detail in Samadi et al. ( 20081 ). 



the computed mode surface velocities, v s , depend significantly on the choice of 
the height, h, in the atmosphere where the mode masses are evaluated. This is 
illustrated in Fig. [2] f or the case of the star a Cen A. 

A recent work bv iKieldsen et al" ( 2008) allows us to estimate the value for 



an effective r^. Indeed, the authors have found that solar modes measured with 
the UCLES spectrograph have amplitudes slightly smaller than those measured 
by the BiSON network. The instruments of the BiSON network use the potas- 
sium (K) resonance l ine, which is formed at an optical depth ~ 0.013. 
iKjeldsen et al. ( 2008f Ts results then suggest that acoustic modes measured by 



UCLES or an equivalent spectrograph (e.g. HARPS) are measured at an effec- 
tive height (rh) slightly below the formation depth of the K line, i.e. at optical 
depth slightly above t 50 o nm — 0.013. Accordingly, in the case of stellar seismic 
observations we will evaluated the mode masses at that optical depth. A more 
rigorous approach would be to compute an effective mode mass by weighting 
appropriately the different mode masses associated with the different spectral 
lines that contribute to the seismic measure. In order to infer accurate mode 
excitation rates from the stellar seismic data, the mode masses representative 
of the observation technique and the spectral lines of the observed star must be 
derived. However, this calls for further studies. 
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Fig. 1. Left: Solar p-mode excitation rates, V, as a function of frequency and de- 
rived from different instruments. The filled circles correspond to seismic data from 
SOHO/GOLF, the diamonds to seismic data from the BiSon network, and the trian- 
gles to seismic data from the GONG network. Right: Same as left panel. V derived 
from GOLF data multiplied by a factor in order that they match at low frequency the 
V derived from GONG or BiSON. 
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Fig. 2. Mode mass evaluated for the case of a Cen A at different heights h above the 
photosphere. The upper curve corresponds to the photosphere (h — 0) and the lower 
curve to the top of the atmosphere (h — 1000 km). The step in h is 200 km. 



4 Theoretical models 

4.1 The inhomogeneous wave equation 

Most of the theoretical models of stochastic excitation adopt GK's approach. 
It consists to solve first, with appropriate boundary conditions, the equation 
that governs the adiabatic wave propagation (also called the homogeneous wave 
equation). This provides the well-known adiabatic displacement eigenvectors 
(£(r, i)). Then, we include in the wave equation of propagation turbulent sources 
of driving as well as a linear damping. The complete equation (so-called inho- 
mogeneous wave equation) is then solved and the solution corresponds to the 
forced mode displacement, Sr asc (r,t) (or equivalently the oscillation mode ve- 
locity v osc = dSr osc /dt). 

A detailed deriv ation of the solution c an be found in I Samadi fc Goupil ( 2001 . 



SG herefater) or in IChaplin et ail (|2005l) . We recall below the main steps 



Equilibrium quantities are represented with a subscript 0. Each variable /, 
except for the velocity v, is written as the sum of the equilibrium quantity , 
/o and an Eulerian fluctuation, /i, f = fo + fi and we retain terms which are 
linear and quadratic in the variables Pi and p\ and neglect, g\, the gravita- 
tional perturbation Q. Accordingly, one obtains for the perturbed momentum 



Neglecting the perturbation of the gravity corresponds to ICowlingl ()l94lh 's approx- 
imation. This approximation remains valid for modes with a high n radial order. 
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and continuity equations: 
dpv 



f)j , V : (pvv) + VPi - Pl g = (17) 

^ + V.(pv) = 0. (18) 

where w osc is the mode frequency, P, p, v and g denote respectively the gas 
pressure, density, velocity and gravity. 

The perturbed equation of state to second order in a Eulerian description is 
given by: 

Pi = c 2 s pi + a s si + u pp p\ + a ss s\ + a ps pisi (19) 

where s is the entropy, a s — [dP/ds) p , c s = A Po/po denotes the average sound 
speed, A = {d\aP / d\n p) s is the adiabatic exponent and a pp , a ss and a ps are 
the second partial derivatives of P versus s and p. Note that Eq. p^|) assumes 
a constant chemical composition (this is indeed the case in the outer convective 
layers) but also constant ionisation rates. 

The velocity field v is split into a component due to the pulsational displace- 
ment Sr osc and a turbulent component u as 

v = v osc + ii (20) 

Linearisation of Eq. (|17til9|) yields for the velocity field, in the absence of 
turbulence (u — 0), the homogeneous wave equation 

d 2 \ 

q£-L)v osc = (21) 

where L is the linear wav e operator (see its expression in SG). With appropri- 
ate boundary conditions (jUnno et alJ Il989h one recovers the usual eigenvalue 
problem : 

L(£(r,t)) = - Lu 2 osc ti(r,t) (22) 

where ui osc is the mode eigenfrequency and £(r, t) = e -4 " 030 * £(r) is the adiabatic 
displacement eigenvector. 

In the presence of turbulence, the pulsational displacement (Sr osc ) is written 
in terms of the above adiabatic solution £(r, t) and an instantaneous amplitude 
A(t) according to Eq. ([B]). Under the assumption of a slowly varying intrinsic 
amplitude A(t), the velocity (v osc ) is related to A(t) and 5r osc according to 
Eq. ©. 

Differentiating Eq. (Tl7|) with respect to t, subtracting the time averaged equa- 
tion of motion, neglecting non-linear terms in v osc , assuming an incompressible 
turbulence (V.tt = 0) and using Eqs. (fT8|) and (fT9|) yields the inhomogeneous 
wave equation 

/ d 2 \ d 
p Q l — -L \ [v osc ] + T> [v^] = —S-C (23) 
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with 



S =S R + S S (24) 
S R = V : (p uu) - V : {(p uu)) (25) 
S s =-V(a s s t ) (26) 



where St is the Eulerian turbulent entropy fluctuations and a s = (dP/dp) s . 
The terms Sr (Eq. l[2"5"])) and 5s (Eq. (|2"6")l ) are two driving sources, namely the 
Reynolds stress tensor and a term that involves the Eulerian entropy fluctua- 
tions. The last term C in the RHS of Eq. (fT8|) gathers terms that involve p\ as 
well as the second order terms of Eq. (fT9|) . C can in principle contribute to the 
driving. However, one can show that its contribution is negligible compared to 
Sr and Ss (see SG, GK). 

The operator T> in the LHS of Eq. (fTSf involves both the turbulent velocity 
field (u) and the pulsational velocity. This term contributes to the dynamical 
linear damping. 

As we will see later, it is more convenient to decompose the Eulerian entropy 
fluctuations in terms of the Lagrangian ones, that is as: 

where sq is the mean entropy. Accordingly, Ss is such that: 
dS s _ f d 



01 = ~ V [j t (^Ss t )-a s u.V s t j (28) 

where we have dropped the term u .Vso since it does not contribute to the 
driving (GK, see also SG). Integration of Eq. (|28[) with respect to time then 
gives 5 S • 



4.2 General solution 

Substituting Eq. (TTJ) into Eq. (|23). yields, with the help of Eq. (^1]). a differen- 
tial equation for A(i). This latter equation is straightforwardly solved and one 
obtains the solution for A: 

A(t) = -^- [ dt' f d 3 xe^ +l ^ t ' e(x).S(x,t') (29) 
2io osc I J -oo Jv 

where I is the mode inertia (which expression is given in Eq. ([9])) and the spatial 
integration is performed over the stellar volume, V. As the sources are random, 
A can only be calculated in square average, (|A| 2 ). This statistical average is 
performed over a large set of realizations. From Eq. (|2"9"|) and with the help of 
some simplifications as detailed in SG, one finds: 
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with 

p p + co 

C 2 = d 3 x d 3 rdTe- luJ °^ T {£* .S 1 £.S 2 ) (31) 

J V J-oo 

where r\ is the mode damping rate (which can be derived from seismic data), 
/ the mode inertia (Eq. ([9])), Xq the position in the star where the stochastic 
excitation is integrated, V is the volume of the convective region, S represents 
the different driving terms, r, and r are the spatial correlation and temporal 
correlation lengths associated with the local turbulence, subscripts 1 and 2 refer 
to quantities that are evaluated at the spatial and temporal positions [xq — — J] 
and [xq + §, §] respectively, and finally (.) refers to a statistical average. 

According to Eqs. (J!}, ((8|) and (|30|) . the theoretical mode excitation rate, V, 
is then given by the expression: 

V = g (32) 



4.3 Driving sources 



The Reynolds stress tensor (Eq. ([23)) ) was identified early on by Lighthilll ( 19521 ) 
as a source of acoustic noise and then as a source of mode excitation (GK). 
This term represents a mechanical source of driving and is considered by most 
of th e theoretical formulations as the dominant contribution to the mode excita- 
tion dGoldreich fc Kee evl 1977b : Dolginov fc Muslimovl 119841: lBalmforthlll992a : 
IStein fc Nordlundll200lHSamadi et al.ll2003aHChaplin et al.ll2005l ). However, as 



pointed-out by lOsakil ( 1990f ). the first calculations by GK's significantly under- 
estimate the power going to the solar modes compared to the observations. 
In order to explain the mode excitation rates derived from the observations, 



Goldreich et al.1 (|1994 GMK hereafter) identified the Lagrangian entropy fluc- 
tuations, i.e. the term Ss t in Eq. (|28p . as an additional driving source. These au- 
thors claimed that this term is the dominant source of driving. However, GMK 
assumed that entropy fluctuations (s t ) behave as a passive sc alar. A passiv e 
scalar / is a quantity that obeys an equation of diffusion (see e.g. lLesieur 1997 ): 



df df „ , 



XV 2 / 



(33) 



where \ is a diffusion coefficient. As shown by SG, assuming as GMK that 8st is 
a passive scalar leads to a vanishing contribution. On the other hand, SG have 
shown that the term a s u . V s t in the RHS of Eq. (|28|) contributes effectively to 
the mode driving. In SG formulation, the so-called entropy source term is then : 



d_ 

di 



S S = V (a s u.Vs t ) 



(34) 



The term m.Vsj in the RHS of Eq. (|34|) is an advective term. Since it involves 
the entropy fluctuations it can be considered as a thermal source of driving. 
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The source term of Eq. was also identified by GK, but was considered 
as negligib l e. It m ust also be pointed out that the theoretical formalisms by 



Balmforthl (11992a!) and Chapli n et al.l ([2005) did not consider this source term. 
According to Samadi et al.l ( 2003al ). this term is not negligible (about ~ 15 % of 
the total power) but nevertheless small compared to the Reynolds stress source 
term (Sr) in the case of the Sun. 

Finally, as seen in Eq. (f3Tj) . Su and Ss lead to cross terms. However, assuming 
as GMK that St behaves as a passive scalar and an incompressible turbulence 
(i.e. V.u = 0), SG have shown that the crossing term between S-ji and Ss 
vanishes. Hence, in the framework of those assumptions, there is no canceling 
between the two contributions (but see Sect. [5]). 



4.4 Length scale separation 

As seen in the RHS of Eq. (l3Tj) . the eigen-displacement £(r) is coupled spatially 
with the source function, S. In order to derive a theoretical formulation that 
can be evaluated, it is necessary to spatially decouple £(r) from S. This is 
the reason why all theoretical formulations explicitly or implicitly assume that 
eddies that effectively contribute to the driving have a characteristic length scale 
smaller than the mode wavelength. Indeed, provided this is the case, £(r) can be 
removed from the integral over r and r that appears in the RHS of Eq. ([30]) (see 
SG). This assumption is justified for low turbulent Mach numbers M t (M t oc 
u/c s where c s is the sound speed). However, at the top of the solar convective 
zone, that is in the super-adiabatic region, M t is no longer small (M t ~ 0.3). 
Furthermore, for G and F stars lying on the main sequence, M t is expected to 
incr ease with the effectiv e temperature and to reach a maximum for M ~ 1.6 Mq 
(see Houdek et al. 19991 ). Hence, for F type stars, significantly hotter than the 



Sun, the length scale separation becomes a more questionable approximation 
(see the discussion in Sect. ITT)) . 



4.5 Closure models 

The second integral in RHS of Eq. (f5U|) involves the term (Si S2) , which is a two- 
point spatial and temporal correlation products of the source terms. Hence, the 
Reynolds stress source term (Eq. ([25])) leads to the two-point correlation product 
of the form ((uu)i («n)j), In the same way, the entropy source term (Eq. (l34|) ) 
leads to the two-point correlation product of the form {(us t )i (ust^)- In both 
case, we deal with fourth-order two-point correlation product involving turbulent 
quantities (that is u and St). Fourth-order moments are solutions of equations 
involving fifth-order moments. In turn, fifth-order moments are expressed in term 
of six-order moments ... and so on. This is the well known closure problem. A 
simple closure model is the quasi-normal approximation (QNA hereafter) that 
permits one t o express four th order moments in term of second order ones (see 
details in e.g. lLesieuilll997l ). that is : 



((uiiij)i {u k ui) 2 ){r,T) = ((ujUj)i) {{u k ui) 2 ) + {{ui)i (1^)2) (( u j)i (ukh) 

+ (« fc )2> <(«i)i («l)a> (35) 
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The decomposition of Eq. (|35p is strictly valid when the velocity is normally 
distributed. The first t erm in the RHS of Eq. (f35|) cancels the term (uu) in 
Eq. (|25l) (see details in Chaplin et alJl2005h . An expression similar to Eq. (j35j) 



is derived for the correlation product ((us t )i {us t )2) (see SG). 
4.6 Adopted model of turbulence 

It is usually more convenient to express Eq. (|35[) i n the frequency (u>) and 
wavenumber (k) domains. We then define tfrij as the temporal and spatial Fourier 
transform of ((14)1 (itj^). For an inhomogeneous, incompressible, isotropic and 
stationary turbulence, there is a relat ion between <f)ij and the kinetic energy 
spectrum E, which is l|Batchelorlll970h : 



where k and ui are the wavenumber and frequency respectively as sociat e d wit h 
the turbulent elements, and 6ij is the Kronecker symbol. Following [Steirj (|1967l ). 
it is possible to split for each layers E(k,u>) as: 

E{k,u)=E(k)xk(fi>) (37) 

where E(k) is the time averaged kinetic energy spectrum and Xk(w) is the fre- 
quency component of E(k,u>). In other words, Xfc( w ) measures - in the fre- 
quency and k wavenumber domains - the temporal correlation between eddies. 
As discussed in Sect. 15.21 the way the eddy time-correlation is modeled has an 
important consequence on the efficiency of the mode driving. A decomposition 
similar to that of Eq. (|3"T|) is performed for the spectrum associated with the 
entropy fluctuations (E s (k,u>)). 

Note that Xk{w) and E(k) satisfy by definition the following normalisation 
conditions: 

/>+oo 

dwX*H = l, (38) 

o 

1 <& 1 

dkE(k) = -( U i) = -(ul)^-ui, (39) 

where u z is the verti cal component of the velocity, = (u 2 ) / (u 2 z ) is the anisotropy 



wnere u z is tne vertical component or tnc velocity, <p = [u )/ \u z ) is tne anisotropy 
factor introduced bv lr*oughl(jl977ft . and Mo is a characteristic velocity introduced 



for convenience. A normalisation condition similar to Eq. (|39[) is introduced for 
E s (k) (see details in SG). 

4.7 Complete formulation 

On the basis of the different assumptions mentioned above, SG then derive for 
radial modes the following theoretical expression for V : 



= — {Cl + CQ (40) 
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where C\ and C| are the turbulent Reynolds stress and entropy contributions 
respectively. There expressions are (see SG): 



Cl 



c% 



4tT 3 £ 
4tT 3 H 



M 



dm po 



S R (m,uj 



M a 2 



dm — ffr(Cr, m) 5s(m, w osc ) 
Po 



(41) 
(42) 



with S# and Ss are the source terms associated with the Reynolds stress and 
entropy fluctuations respectively: 



Sb = 
Ss = 



dk 
dk 



E 2 (k,m) 
E s (k, m)E(k, m 



dwXk{u i 

+ 00 



-w,m)xfe(w,m) (43) 
dwXfc(wo S c + w,mk)i(w,m) (44) 



In Eq. ([41]) and ([42]) , po is the mean density, Q and 7i are two anisotropic factors 
(see their expressions in SG), and finally .g r (£r, to) is a function that involves the 
first and the second derivatives of £ r , its expression is: 



7r(6)TO) 



1 da s d^ r d 2 £ r 
a, G?r dr c?r 2 



(45) 



It is in general more convenient to rewrite Eqs. (j4lj) and (|42|) in the following 
forms: 



c 2 R = A^g 



M 



dm 



Po ut 



k : '> 



0^0 



r 2 - 



dm 



(q s su ) 2 
po w o 



<§ij(m, L0 OSC ) , 



g r (£ r ,m) S s (m,u osc ) 



(46) 
(47) 



where we have defined the dimensionless source functions = (fcg wo / Mq) Sr 
and = (fcp ojo / (uq s)) 5/?,, s and where s is the rms of the entropy fluctuations. 
We have introduced for convenience the characteristic frequency ujq and the 
characteristic wavenumber k ; they are defined as: 



uj = k u 
~A 



ko 



(48) 
(49) 



where A is a characteristic size derived from E{k) and uo is the characteris- 
tic velocity given by Eq. (j39|) . For future use, it is also convenient to define a 
characteristic time To as: 

2tt A 

to = T— = — (50) 



fc uo 



uq 
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From Eq. (|46p we can show that the driving by the Reynolds stress is locally 
proportional to the kinetic energy flux. Indeed, the flux of kinetic energy in the 
vertical direction is by definition: 



F kin = w S kin = w [ - po u z ) = - \j — p ul , 



3 /3 

2 Vi' 



(51) 



where -Ekin = (1/2) PqU 2 is the kinetic energy per unit volume. Substituting 
Eq. ([51]) into Eq. j46]) yields the relation: 



Clcx 



dmFkin A 4 



dr 



S R (m,u osc ) 



(52) 



Concerning the driving by the entropy fluctuations, we can show that locally 
this driving does not only depend on Fkm but also on the convective flux (F c ). 
Indeed, lets define as GMK the quantity: 



K = 



Substituting Eq. ([53]) into Eq. ([47]) yields the relation: 



Cloc 



dmFkin A 4 V?T 2 



Wosc 



Ss(m,u; osc ) , 



(53) 



(54) 



where we ha ve defined as in SG the quantity T 2 = A 2 g T . Finally, since 1Z oc 
F c /F^ a (see Samadi et al. 20061) . we can conclude that locally the driving by 
entropy source term is proportional to Fkin and to the square of the ratio 1Z oc 
F c /Fki n . 



5 Turbulent spectrum 

As seen in Sect. 14.6] the model of stochastic excitation developed by SG involves 
E(k,u), the turbulent kinetic spectrum as well as E s (k,oj), the spectrum associ- 
ated with the turbulent entropy fluctuations. Both spectra are split in terms of a 
time averaged spectrum (E(k) for the velocity and E s (k) for the entropy fluctu- 
ations), and a frequency component Xfc( w ) ( see Sect. 14. 6[ ). Different prescriptions 
were investigated for both components. The results of these investigations are 
summarized in Sect. 15. II for E(k) and in Sect. 15.21 for Xfe(w). 



5.1 Time averaged spectrum, E(k) 

Two approaches are commonly adopted for prescribing E{k). The classic one 
is to assume an analytical function derived either from theoretical considera- 
tions or empirical ones. The more commonly use d analytical spectrum is the so- 
called Kolmogorov spectrum (|Kolmogorovlll94lh . which derives originally from 
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Oboukhovl (|1941r )'s postulate that energy is transferred from the large scales 
to the small scales at a con stant rate. O ther theoretical spectra, such as the 



so-called Spiegel's spectrum (Spiegel! 1962 ). or purely empirical spectra, such as 
those proposed bv lMusielak et aL (1994), were also considered. All of these an- 



alytical functions differ from each other by the way E(k) varies with k. But for 
all of them, it is required to set a priori the characteristic wavenumber, kg, at 
which energy is injected into the turbulent cascade. The second approach con- 
sists to obtain E(k) directly from hydrodynamical 3D simulations. This method 
has two advantages: it provides both the k dependence of E{k) and the charac- 
teristic wavenumber k$. On the other hand, the inconvenient is that such method 
depends on the quality of the 3D hydrodynamic al simulation. 



These two approaches have been compared in lSamadi et al.l (|2003bl ). Among 
the different analytical functions tested, the best agreement with a solar 3D 
simulati on was found wi t h the so-called "Extended Kolmogorov Spectrum" de- 
fined bv iMusielak et aD ( 1994 ). This spectrum increases at low scales as k +1 
and decreases at low scales according to the Kolmogorov spectum, i.e. as fc -5 / 3 . 
However, due to the limited spatial resolution of the solar simulation used, the 
Kolmogorov scaling is validated over a limited range only. Nevertheless, the ma- 
jor part of the excitation arises from the most-energetic eddies, also refered to the 
energy bearing eddies. Accordingly, the contribution of the small scales, that are 
not resolved by the present 3D simulations, are expected to be relatively small. 
However, to confirm this, a quantitative estimate must be undertaken. 

More important is the choice for the characteristic wavenumber fco- Indeed, 
the integrands of Eq. ([46]) and Eq. (|47|) are both proportional to fc^ 4 . Accord- 
ingly, the computed V are very sensitive to the choice for ko. This characteristic 
wavenumber can be obtained from 3D simulations. However, by default, one 
usually relates fco to the mixing-length /Imlt according to: 



fc 



• MLT 



2tt 



0A 



MLT 



(55) 



where vImlt = ex H p is the mixing- length, a the mixing-length parameter, H p 
the pressure scale height, and (3 a free parameter, whi ch is usually s et to a 
value of the order of one. The solar 3D simulation used by Samadi et al. ( 2003bh 
indicates that in the Sun k ~ 3.6 Mm" 1 at the top of the excitation region. This 
characteristic wavenumber corresponds to an horizontal size of the granules of 
A g = 2ir/ko ~ 2 Mm. This horizontal size is reached at the top of the excitation 
region with a value of (3 that depends on the adopted value for a and the solar 
ID model used. For other stars, 3D simulations are rarely available. In that 
case, one usually assumes for (3 the same value that the one adopted for the Sun. 
Hence, an open and important question is whether or not the parameter /3 can 
be kept the same for other stars as for the Sun. 



5.2 Eddy time-correlation, Xk(<+>) 

Most of the theoretica l formulations explicitly or implicitly assume a Gaus- 
sian function for Xfe(w) (jGoldreich fc Keelevlll977bl : lDolginov fc Muslimovill984 : 
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Gold reich et al ]ll994t lBalmfortblll992al ISamadi et afkOQll IChaplin et ap2005l) . 



However, 3D hydrodynamical simulations of the outer layers of the Sun show 
that, at th e length associated with the energy bearing eddies, \k is rather 
Lorent zian (Samadi et al.l 2003al ). This is well illustrated in Fig. [3j As pointed- 
out by IChaolin et alJ (j2005h . a Lorentzian Xk is also a result predicted for the 
largest, mo st-energ e tic ed dies by the time-dependent mixing- length formulation 
derived by iGough (1977). Therefore, there is some numerical and theoretical 
evidences that Xk is rather Lorentzian at the length scale of the energy bearing 
eddies. 

As shown bv lSamadi et al. ( 2003a[ ). calculation of the mode excitation rates 
based on a Gaussian Xk results for the Sun in a significant under-estimation of 
the maximum of V whereas a better agreement with the observa tions is found 
when a Lorentzian Xk is used. A similar conclusion is reached bv lSamadi et all 
(2008) in the case of the star a Cen A. These results are illustrated in Fig. [5] in 
the case of the Sun and in Fig. [6] in the case of a Cen A. 




Fig. 3. Eddy time-correlation function, Xki as a function of frequency v for the layer 
where the radial component of the velocity is maximum. The filled dots represent 
~ Xk obtained from a solar 3D simulation with an horizontal resolution of ~ 25 km 
(Sa madi et alJ l2003a[ ) . Xk is shown here for the wavenumber fc at which E(k) peaks. 
The solid line represents a Lorentzian function and the dashed line a Gaussian function. 



The excitation of low- frequency modes (y < 3 mHz) is mainly due to the large 
scale eddies. However, the higher the frequency the more important the contri- 
bution of the small scales. 3D solar simulations show that, at small scales, Xk 



is neither Lorentzian nor G aussian (jGeorgobiani et aL 2006f ). Hence, according 



to Georgobiani et al. ( 20061 ). it is impossible to separate the spatial component 
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Fig. 4. Top : Same as Fig. [3] for a solar 3D simulation with an horizontal resolution 
of ~ 50 km l|Belkacem et al.ll2006al ). Bottom: Same as top for a solar 3D simulation 
with an horizontal resolution of ~ 120 km. 
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E(k) from the temporal component at all scales with the same simple analyti- 
cal functions. However, such results are obtained using Large Eddy Simulation 
(LES). The way th e small scales are treated in LES can affects our description of 
turbulence. Indeed, He et a l. (2002) have shown that LES results in a Xfc( w ) that 
decreases at all resolved scales too rapidly with w with re spect to direct numer- 
ical simulations (DNS). Moreover. I Jacoutot et al. (2008b) found that computed 
mode exc itation rates d e pend significantly on the adopted sub-grid model. Fur- 
thermore, Samadi et al.l ( 2007 ) have shown that, at a given length scale, \k tends 
toward a Gaussian when the spatial resolution is decreased. This is illustrated in 
Fig. [4] by comparison with Fig. [3l In summary, the numerical resolution or the 
sub-grid model can substantially affect our description of the small scales. Im- 
proving the modeling of the excitation of the high frequency modes then requires 
more realistic and more resolved hydrodynamical 3D simulations. 

Up to now, only analytical functions were assumed for Xk(uj). We have here 
implemented, for the calculation of V, the eddy time-correlation function derived 
directly from long time series of 3D simulation realizations with an intermediate 
horizontal resolution (~ 50 km). As shown in Figs. [5] and [SI the mode excita- 
tion rates, V, obtained from x| D , are found to be comparable to that obtained 
assuming a Lorentzian \k, except at high frequency in the case of the Sun. This 
is obviously the direct consequence of the fact that a Lorentzian \k reproduces 
rather well xjP (see Fig. [3]), except at high frequency where Xfc D decreases more 
rapidly than the Lorentzian function (see Fig. 2] left). At high frequency, cal- 
culations based on a Lorentzian \k result in larger V and reproduce better the 
helioseismic constraints than those based on xjf ( see Fig. EJ. This indicates 
perhaps that Xk decreases more rapidly with frequency than it should. This 
is consistent with iHe et al.l (|2002l ) 's results who found that LES predict a too 
ra pidly decrease with v compared to the DNS (see above). 

Chaplin et al. I (|2005l) also found that the use of a Gaussian \k severely under- 
estimates the observe d solar mode excitation rates. However, in contrast with 
Samadi et al. ( 2003al ). they mention that a Lorentzian Xfc results in a severe 



ov er-estimation f o r the low-frequency modes. In order to illustrate the results 
bv lChaplin et~aj] (|2005h . we have computed the solar mode excitation rates us- 
ing their formalism and a solar e nvelope equilibrium model similar to the one 
considered by these authors (see Samadi et al. 2003bl ). The result is shown in 
Fig. [Jj We clearly see that the mode excitation rates computed using a Gaus- 
sian Xk overestima te by ~ 20 th e seismi c constraints. This result is consistent 
with this found bv | Samadi et al.l (|2003al ). On the other hand, in contrast with 
Samadi et al.l (|2003al ). the modes with frequency below v ~2 mHz are severely 



over-estimated when a Lorentzian Xk is assumed. It should be pointed out that 
the excitation of modes with frequency v < 2 mHz occu rs in a region more ex - 
tended than covered by the solar 3D simulation used by [S amadi et al.l (|2003bl ). 
On the other hand, the pure ID modeling by IChaplin et al.1 (|2005l ) , includes 
all of the convective zone. The severe over-estimation at low frequency of the 
mode excitation rates, is explained by the authors by the fact that, at a given 
frequency, a Lorentzian \k decreases too slowly with depth compared to a Gaus- 
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1 2 3 4 5 

v [mHz] 



Fig. 5. Solar p-mode excitation rates as a function of v. Filled circles and diamonds 
correspond as in Fig.[T]to seismic data from SOHO/GOLF and BiSON network respec- 
tively. The lines correspond to semi-theoretical calculations based on different choices 
for Xk'- Lorentzian \u (solid line), Xfc3D i.e. Xfc derived directly from the solar 3D 
simulation (dashed line), and a Gaussian \k (dot-dashed line). 



a C^n A 

Lorentzian 
Gaussian Xk 




Fig. 6. Same as in Fig. [5] for the case of a Cen A. 
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sian Xk- Consequently, for the low-frequency modes, a substantial fraction of the 
integrand of Eq. (|4Tj) arises from large eddies situated deep in the Sun. This 
might suggest that, in the deep layers, the eddies that contribute efficiently have 
rather a Gaussian Xk- However, this remains an open issue. 




1 2 3 4 5 

v [mHz] 



Fig. 7. Same as Fig . [5] The lines correspond to calculations using the formalism by 
IChaplin et all(|2005l ). Two choices for \k was considered : a Lorentzian \k (solid line) 
and a Gaussian Xfc (dashed line). In both calculations, driving due to the entropy 
fluctuations is not included. 



6 Closure models and anisotropy 



The decomposition of Eq. (|35[) assumes the quasi-normal approximation (QNA). 
However, it is well known that the departure from the QNA is important in a 
strongly turbulent medium. In addition, a closure model base d on the QNA 
does not ensure the positiveness of the energy (see details in e.g. lLesieurlll9971 ). 
Furthermore, the QNA is strictly valid only for normally distributed fluctuat- 
ing quantities with zero mean. However, the upper-most part of the convection 
zone is a turbulent convective medium composed of essentially two flows that are 
asymmetric with respect to each other. Hence, in such a medium, the probability 
distribution function of the fluctuations of t he vertical veloci t y and t emperature 
do not foll ow a Gaussian law. As ver ified by iBelkacem et al.1 ( 2006al B06a here- 
after) and lKupka fc Robinson! (|2007l ). departure from the QNA is important in 
the upper part of the solar convective zone. Indeed, this approximation under 
estimates, in the quasi-adiabatic region, by rs 50 % the fourth-order moment of 
the vertical velocity derived from a solar 3D simulation. 
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The term in the LHS of Eq. (|55|) corresponds to a two-point correlation 
product involving the velocity, i.e. ((ui Uj)\ (uk 111)2) (r, r) where r and r are the 
spatial correlation and temporal correlation lengths respectively. For r — > and 
t — > 0, this term reduces to a one-point correlation product, (uiUj u^ui) , also 
refered as a fourth-order moment (FOM hereafter). As we shall see below, it is 
possible to derive an improved closure model for this term that does not rely 
on the QNA. However, we still require a prescription for the two point corre- 
lation products involving the velocity ({(uiUj)i (u^ 111)2) (r,r)) and the entropy 
fluctuations ({(u s t )i (u s t )2)(r, r)). For radial modes or low t order modes, only 
the radial component of the velocity (w) matters. Hence, for these modes we 
require a presc r iption for (wf w 2 )(r, t) and ((w St)i (w St)a)(r, r). By default, 



Belkacem et al.l (|2006b[ ) have proposed that {w\ ra^r, t) varies with r and r in 



the same way than in the QNA (Eq. (f35|) ). that is: 

{w\ wf) = ~y {w\ w|) Q na , (56) 

where IC W is a constant and (w 2 w 2 )qna is the two-point correlation product 
given for w according to the QNA (Eq. ([35)0 . Accordingly, the contribution of 
the Reynolds stress (C R} Eq. (|4T|) ) is modified as: 

C 2 R = 4ir 3 g J dm Po (^j ^S R (m,uj osc ) (57) 

Note that the contribution of the entropy fluctuations (C|, Eq. (|42]0 still assumes 
the QNA. This inconsistency has a small impact on computed mode excitation 
rates since C| is significantly smaller than C R , at least for stars that are not too 
hot (but see Sect. El). 

The constant IC W is determined in the limit case where r — > and r — > 0. 
Indeed, when r — » and r — » 0, we have: 

(w4) = ^ (w 4 )qna ( (5g) 

where (w 4 ) is by definition the fourth-order moment (FOM hereafter) associated 
with w and (w 4 )qna is the one given by the QNA. In the same way, Eq. (I35|) 
gives : 

(w 4 )qna = 3 (w 2 ) 2 . (59) 
Using Eqs. ([58)) and ([59)) . we then derive the constant 1C W : 

(w 4 ) Q na {w z Y 

which is by definition the Kurtosis. This quantity measures the oblateness of the 
probability density function (see e.g. B06a). For normally distributed w we have 
JC W = 3. The Kurtosis then measures the departure of the FOM from the QNA. 
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Closure models more sophi sticated than the QNA can be used. Amo ng those, 

the tw o-scale mass flux model ( Abdella fe McFarlandll997 ) improved bv lGrvanik fe Hartmann 



(2002) takes the asymmetries in the medium into account but is only applica- 



ble for quasi-laminar flows. For K. w , Grvanik fe Hartmann! ( 20021 ) obtained the 
following expression: 

IC w = (l + Sl) (61) 
with the skewness, S w , given by: 

Sw 55 = (62) 

where a is the mean fractional area occupied by the updrafts in the horizontal 
plane. In the QNA limit, i.e. when the random quantities are distributed accord- 
ing to a Normal distribution with zero mean, we necessarily have S w = 0. Hence, 
in the QNA limit, Eq. (IFTj) d oes not match the expected value i.e. JC W = 3. Then, 



Grvanik fc Hartmann! (|2002l ) proposed to modify Eq. (|6ip as follows: 

1C W =3(1 + ±S 2 J. (63) 



Figure M shows that the FOM based on Eq. (JB3J with S w given by Eq. ([62]). 
results in a negligible improvement with respect to the QNA. However, when S w 
is derived directly from the 3D simulation and plugged into Eq. (|63|) . Eq. ([63]) 
is then a very good evaluation of the FOM de rived from a 3D simuation of the 
ou ter layer of the Sun as v erified by B06a a ndlKupka fc Robinson! 



Belkacem et al.l (|2006al ) have generalized iGrvanik fc Hartmann! (|2002l )'s ap- 



proach by taking the skewness introduced by the presence of up- and down-drafts 
and the turbulent properties of each flow into account. Accordingly, they have 
derived a more accurate expression for S w (see the expression in B06a). As shown 
in Fig. [H calculations of the FOM based on Eq. (|63|) and their expression for S w 
reproduce rather well - in the quasi-adiabatic region - the FOM derived from 
t he solar 3D simulation. 

iBelkacem et~aH (|2006bl) have computed mode excitation rates, V according 



to Eq. (|57|) with the Kurtosis 1C W given by Eq. (|63[) and with the skewness 
S w computed according to B06a's closure model. The maximum in V is found 
about 30% larger than in calculations based on the QNA and fits better the 
maximum in V derived from the helioseismic data. This increase is significantly 
larger than the entropy contribution (the term Ss in Eq. (|42p. which is of the 
order of ~ 15 %, see Sect. [8]). We stress that, however, 30 % is of the same order 
as the difference between seismic constraints of different origins (SOHO/GOLF, 
GONG, BiSON). These results are illustrated in Fig. [9] 



24 R. Samadi 



2.0 




Fig. 8. Fourth-order moment (FOM) of the velocity, (w4) = 1C W (u>2) 2 , as a function 
of depth z, normalized to the FOM derived from the 3D simulation. In all cases the 
Kurtosis K, m , (Eq. J60J) is calculated according to Eq. (163[) but with different skewness, 
Sw The solid line S w is compu t ed ac cording to B06a's closure model, the dashed line 
assumes iGrvanik fc Hartmannl (|2002l Vs expression for S w (Eq. ()62p ) and finally the 
dotted line assumes the QNA, that is Sw — and K, 



7 Importance of the stellar stratification and chemical 
composition 

7.1 Role of the turbulent pressure 

Rosenthal et~aTl (|l999f ) have shown that taking the turbulent pressure into ac- 
count in a realistic way in the ID global solar models results in a much better 
agree ment between observed and theoretical mode frequencies of the Sun. Fol- 
lowing Rosenthal et al. ( 19991 ). Samadi et al. ( 2008 ) have studied the importance 
for the calculation of the mode excitation rates of taking the turbulent pressure 
into account in the averaged ID model. For this purpose, they have built two 
ID models representative of the star a Cen A. One model (here refered as the 
"patched" model), has its surface layers taken directly from a fully compressible 
3D hydrodynamical numerical model. A second model (here refered as "standard" 
model), has its surface layers c omputed using standa rd physics, in particular con- 
vection is described according lBohm-Vitens3(ll958l) 's mixing-length local theory 
of convection (M L T) an d turbulent pressure is ignored. 

Sa madi et al. (2008) found that the calculations of V involving eigenfunc- 
tions computed on the basis of the "patched" global ID model reproduce much 
better the seismic data derived for a Cen A than calculations based on the eigen- 
functions computed with the "standard" stellar model, i.e. built with the MLT 
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Fig. 9. Same as in Fig. [5] The thick lines correspond to calculations where the Reynolds 
stress contribution is computed according to Eq. (|57|l . The Kurtosis (K w ) is computed 
here in a different manner: for the solid line K, w is obtained directly from a 3D solar 
simulation, for the dashed line the Kurtosis is calculated according to Eq. (|63[) where 
the skewness (S w ) is obtained from B06a's closure model, and finally for the dot-dashed 
line we have assumed the QNA, that is S w = and K, w = 3. 



and ignoring turbulent pressure. This is because a model that includes turbulent 
pressure results in lower mode masses M. than a model that ignores turbulent 
pressure. This can be understood as follows: Within the super-adiabatic region, 
a model that includes turbulent pressure provides an additional support against 
gravity, hence has a lower gas pre ssure and density than a model that does not in - 
clude turbulent pressure (see also lNordlund fc Steinlll999l : Rosenthal et al.llT999l ) . 



As a consequence, mode inertia (Eq. ([9|)) or equivalently mode masses (Eq. (fTTj)) 
are then lower in a model that includes turbulent pressure. 

7.2 Role of the surface metal abundance 

Samadi et al. ( 2009(J ) have recently studied the role of the surface metal abun- 



dance on the efficiency of the stochastic driving. For this purpose, they have 
computed two 3D hydrodynamical simulations representative - in effective 
temperature and gravity - of the surface layers of HD 49933, a star which 
is rather metal poor compared to the Sun since its surface iron-to-hydrogen 
abundance is [Fe/H]=-0.37. One 3D simulation (hereafter labeled as SO) has a 
solar metal abundance and the other (hereafter labeled as SI) has [Fe/H] ten 
times smal ler. For each 3D simu lation they have build a "patched" model in the 



manner of Samadi et al. (2008) and computed the acoustic modes associated 



with the "patched" model. 
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As seen in Fig.flUJ the mode excitation rates V associated with SI are found 
to be about three times smaller than those associated with SO. This difference 
is related to the fact that a lower surface metallicity results in a lower opacity, 
and accordingly in an higher surface density. In turn, the higher the density, the 
smaller are the convective velocities to transport by convection the same amount 
of energy. Finally, smaller conve ctive velocities result in a less efficient driving (for 



details sec Samadi et al. 2009c). This conclusion is qualitatively consistent with 



that bv lHoudek et all (|l999) who - on the basis of a mixing-length approach - 
also found that the mode amplitudes decrease with decreasing metal abundance. 




Fig. 10. Mode excitation rates, V, as a function of the mode frequency (y ) obtained for 
two 3D models with the effective temperature and the surface gravity of HD 49933 b ut 
with two different surface metal abundances (see Sect. I7.2l and Samadi et al.l (|2009d )). 
The solid line corresponds to the 3D model with the metal abundance (SO) and the 
dashed line to metal poor 3D model (SI). The dot-dashed line correspon ds to the mode 
excitat ion rates derived for the specific case of HD 49933 as explained in iSamadi et al.l 
l|2009d ). 



Using the s eismic determinations of the mode linewidths measured by CoRoT 
for HD 49933 (jBenomar et all [2009) and the theoretical mode excitation rates 
computed for the specific case of HD 49933. ISamadi et all (|2009bf) have derived 
the theoretical mode amplitudes of the acoustic modes of HD 49933. Except at 
rather high frequency (y > 1.9 mHz), their amplitude calculations are within 
approximately 1-a in agreement with the mode amplit udes derived from the 
CoRoT data (for more details see ISamadi et al.l l2009bh . They also show that 
assuming a solar metal abundance rather than the observed metal abundance of 
the star would result in larger mode amplitudes and hence in a larger discrepancy 
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with the seismic data. This illustrates the importance of taking the surface metal 
abundance of the solar-like pulsators into account when modeling the mode 
excitation. 



8 Contribution of the entropy fluctuations 



Using the method summarised in Sect. 19.21 IStein fc Nordlund ( 2001 ) have com- 
puted directly from a 3D simulation of the surface of the Sun the contribution of 
the incoherent entropy fluctuations (Eq. ([26]) V They also found that the entropy 
fluctuatio n is small compared to the Reynolds stress contribution. However, as 
shown by ISamadi et alJ <|2007t ) , the relative contribution of the entropy to the 
total excitation rate increases rapidly with the effective temperature, T e g. For 
instance, the solar-like pulsat or HD 49333 has a significantly higher T e g than 
the Sun. ISamadi et al.l (|2009d l found that for this star the entropy fluctuations 
contributes up to ~ 30 % while it is only about 15 % in the case of the S un (see 
Samadi et al.ll2007l ) a nd in the case o f a Cen A fsee lSamadi et al.l 120081 ). 

As poi nted-out bylHoudekl (|2006l ). the solar and stellar 3D simulations per- 
formed bv lStein et al.l (j2004l ) show some partial canceling between the Reynolds 
stress contribution (Sr, Eq. (|25"|) ) and contribution due to the entropy (Ss , 
Eq. ([26]) ). This cancellation increases with increasing T e g (see Stein et all 2004 ). 
In the theoretical model of stochastic excitation, the cross terms between the 
entropy fluctuation s and th e Reyn olds stresses vani sh (see Sect. 14.31). As origi- 
nally suggested by Houdek ( 2006t ) and discussed in Samadi et al. ntioOQbJ) . the 
existence of a partial canceling can decrease the mode amplitude and improve 
the agreement with the seismic observations. However, there is currently no the- 
oretical modeling of the interference bet ween these two terms (see the discussion 
in Sect, fill and in lSamadi et al.1 (|2009fJ ). 



9 Alternative approaches 
9.1 Energy equipartition 

Under certain conditions that we will emphasize below, GK have shown that 
there is an equipartition of kinetic energy between an acoustic mode and the 
resonant eddy. To derive this principle, GK assume that the acoustic modes 
are damped by turbulent viscosity and excited by the Reynolds stresses. We 
reproduce here their demonstration. For the sake of simplicity, we will consider 
modes with w osc To < 1 where To is the characteristic time of the energy bearing 
eddies typically located in the upper part of the convective zone, that is the 
region where the driving is the most vigorous. Furthermore, we neglect as did 
GK the driving by the entropy fluctuations (Eq. (|42[) ). According to Eqs. (1401) 
and (|52p . we have roughly for acoustic modes with uj osc tq < 1: 

2 

E eddy A u , (64) 



V oc 



j dm 


d^ 




dr 
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where A is the characteristic size of the energy bearing eddies, uq their char- 
acteristic velocity (Eq. ([39]) ), r = A/uq their characteristic lifetime (Eq. ([50])). 
and i^cddy = (3/2) po Uq A 3 their total kinetic energy. Let fc sc be the vertical 
oscillation wave number. We have then d£ r /dr — ik osc £r- We further assume 
that - in the driving region - the acoustic waves are purely propagating. This 



assumption then implies ui osc - 
we can simplified Eq. (|64|) as: 



fc osc c s where c s is the sound speed. Accordingly, 



V oc 



OSC 

I 



dm 



E, 



eddy 



Au . 



(65) 



In the region where the mode are excited, -Eeddy, uq, and c s vary quite rapidly. 
However, again for the sake of simplicity we will assume that these quantities are 
constant and evaluate them at the layer where the excitation is the most efficient, 
i.e. at the peak of the super-adiabatic temperature gradient. The integration of 
Eq. ([64]) can be approximated as 



V oc 



1 / u>, 



E eddy Au I dm£, 



(66) 



Using the expression of the mode inertia (Eq. (J9j) ) , we can finally simplify 
Eq. (H| as: 



V oc 



-Eeddy A Uq 



(67) 



Modes damped by turbulent viscosity have their d amping rates rj given by 
( Ledoux fc Walravenlll958HGoldreich fc Keelevlll977al ). 



dr 



(68) 



where v t is the turbulent viscosity. The simplest prescription for v t is the con- 
cept of eddy- viscosity. This consists in assuming v t = w A = r u$. Obviously 
the turbulent medium is characterized by eddies with a large spectrum of size. 
However, only the eddies for which lu osc rj « 1 are expected to efficiently damp 
the mode with frequency w osc . Since we are looking at the modes such that 
^osc t\ ;$ 1, only the largest eddies efficiently damp the mode, that is the eddies 
with size A. Accordingly, we adopt v t = uq A. With the same simplifications and 
assumptions as those used for deriving Eq. (|67[) . we can simplify Eq. (|68p as: 



rj oc 



^osc 
Cs 



Au 



(69) 



From Eqs. ([5]), ff^Tj) and ([59|1, we then derive the mode kinetic energy: 



-^osc -^eddy 



(70) 
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Eq. ([70]) highlights an equ ipartition of kinetic energy between an aco ustic mode 
and the resonant eddies. IChristensen-Dalsgaard fe Frandsenl (|1983l ) used this 
"equipartition principle" to derive the first quantitative estimate of solar-like 
oscillations in stars. The relation of Eq. (I70|) was derived by assumin g that 
modes are damped by turbulent viscosity. However, as pointed-out by lOsaki 
dl990l). theoretical mode li ne-widths, r — 77/71", computed in the manner of 
Goldreich fc Keelevl (|l977eJ) . i.e. ass uming a viscou s damping, are underesti- 



mated compared to the observations. iGouehl (119761 proposed a different pre- 
scription for v t . Nevertheless, assum ing [Goughl ( 1976h 's prescription a l so res ults 
in similar r (see Balmforth 1992b). On the other hand, Xiong et al. ( 2000h re- 
port that the turbulent viscosity is the dom i nant s ource of damping of the radial 
p modes. As discussed recently by iHoudek ( 20081 ). there is currently no consen- 
sus about the physical processes that contribute dominantly to the damping of 
p modes. If the damping due to turbulent viscosity turns out to be negligible, 
then there is no reason that the balance between the mode kinetic energy and 
the kinetic energy of resonant eddies holds in general. 



9.2 "Direct" calculation 



The model of stochastic excitation presented in Sect.|4]is based on several simpli- 
fications and assumptions concerning th e turbulence and the sourc e terms. There 
is an alternative approach proposed by Nordlund fc Steir] ( 2001 ) that does not 
rely on such simplifications and assumptions. In such approach, the rate at which 
energy is stochastically injected into the acoustic modes is obtained directly from 
3D simulations of the outer layers of a star by computing the (incoherent) work 
performed on the acoustic mode by turbulent convection. In their approach, the 
energy input per unit ti me into a given aco u stic m ode is calculated numerically 
according to Eq. (74) of iNordlund fc Steinl (|2001l ) multiplied by S, the area of 
the simulation box, to get the excitation rate (in J s _1 ) : 



"P3d(w OS c) 



.S 



■ Av 



dr AP nad (r, w osc ) 

or 



(71) 



where AP^^r, uj) is the discrete Fourier component of the non-adiabatic pres- 
sure fluctuations, AP n& d(r, t), estimated at the mode eigenfrequency uj osc = 
2ttvq 1 £ r is the radial component of the mode displacement eigenfunction, Av = 
1/T S the frequency resolution corresponding to the total simulation tim e T s and 
£uj n .r. is the normalised mode energy per unit surface area defined in lNordlund fc Stein 



(|200ll their Eq. (63)) as: 



2 

osc 



dr 



*'(s) 



(72) 



Eq. (|71"j) corresponds to the calculation of the PdV work associated with 
the non-adiabatic gas and turbulent pressure (Reynolds stress) fluctuations. In 
contrast to the pure theoretical models (see Sect. 2]), the derivation of Eq. (|71[) 
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does not rely on a simplified model of turbulence. For instance, the relation of 
Eq. (|36[) is no longer required. Furthermore, they do not assume that entropy 
fluctuations behave as a passive scalar (Eq. ([33]) ). However, as for the theoretical 
models, it is assumed that £ r varies on a scale-length larger than the eddies 
that contributes effectively to the driving (this is the so-called "length-scale 
separation", see Sect. 14.41) . In addition, Eq. (|7Tj) implicitly assumes the quasi- 
Normal approximation (Eq. (|35| V 

The expression of E q. ([7T|) has been applied to the case of the Sun by 
Stein fc Nordlundl < 200lh . These authors obtain a rather good agreement be- 
tween V 3 v (Eq. ([TT jO and the solar mode excitation rates derived from the 



GOLF instrument 
tion rates derived 



by 
by. 



Roca Cortes et al 



., .Stein fc Nordlund 
Roca Cortes et al.l (|l999F 



are ■ 



( 1999h . However, solar mode excita- 
(|200Ir ) from the seismic analysis by 
for a reason that remains to be und erstood - sys- 
temat i cally lower than th ose derived from the seismic analysis bv lBaudin et al.l 
(|2005L lstein et all (|2004l ) have computed V 3D (Eq. [7T]) for a set of stars located 
near the main sequence from K to F and a subgiant K IV star. The comparison 
between these calculations an d those based on SG's formalism has been under- 
taken bv lSamadi et~aH (|2007l ). The maximum in V3B was found systematically 
lower than those from calculations based on S G's formalism (Eqs. ([40 |) -(j44 | ). 
These systematic differences were attributed by ISamadi et al. ( 2007T ) to the low 



spatia l resolution of the hydrodynamical 3D simulations computed bv lStein et al.l 
(|2004h . 



10 Stochastic excitation across the HR diagram 
10.1 Mode excitation rates 



Using several 3D simulations of the surface of main sequence stars. ISamadi et al 



(2007) have shown that the maximum of the mode excitation rates, 'Pmax, varies 
with the ratio L/M as (L/M) a where L and M are the luminosity and the mass 
of the star respectively and a is the slope of this scaling law. Furthermore, they 
found that the slope a is rather sensitive to the adopted function for Xk'- cx=3.1 
for a Gaussian Xk and a=2.6 for a Lorentzian one. 

The increase of V max with L/M is not surprising: It should first be noticed 
that, even though the ratio L/M is the ratio of two global stellar quantities, it 
nevertheless essentially characterizes the properties of the stellar surface layers 
where the mode excitation is located since L/M oc T^ s /g. Indeed, by definition 
of the effective temperature, T c g , and the stellar radius R, the total luminosity 
of the star, L, is given by the Steffan's law: L = AiraT^ R 2 where a is Stef- 
fan's constant. Furthermore, the surface gravity is g = GM/R 2 where G is the 
gravitational constant. Accordingly, L/M oc T^ s /g. 

Second, as we will show now, it is possible to roughly explain the dependence 
of Pmax with g and T c 4 ff . Eq. (|67|) can be rewritten as: 



UJr. 



" [^) 1% • (73) 
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where 

3 

*kin = g Po M o ( 74 ) 

is by definition the flux of kinetic energy per unit volumc|l and uq is the charac- 
teristic velocity given by Eq. (|39|) . 

The characteris tic size A is approxi mately proportional to the pressure scale 
height H p (see e.g. ISamadi et al. 20081) . From hydrostatic equilibrium, we have 



P = pg H p . Assuming now the equation of state of a perfect gas, we then derive 
H p oc T I g. The sound speed is given by the relation <? s — F\ P/p. Accordingly, 
using again the perfect gas equation, we then have c| oc T. From these simplifi- 
cations, we can simplify Eq. (|T3[) as: 



V ^uj 2 osc F Mn T 3 g~ A . (75) 

In the framework of the mixing-length approach, it can be shown that Fkin 
is roughly proportional to the convective flux F c . Indeed, in this framework, 
the eddies are accelerated by the buoyancy force over a distance equal to the 
mixing-length A = a H p where a is the mixing-length parameter. Accordingly, 
the kinetic energy o f the eddies, -E ddy, is given by (see the lecture notes by 
Bohm-Vitenselll989| ) 



E cddy = \p ul A 3 = g (ApA 3 ) A (76) 

where Ap is the difference between the density of the eddy and its surroundings. 
In the Boussinesq approximation, the perturbation of the equation of state gives: 

where AT is the difference between the temperature of the eddy and its sur- 
rounding. Now, the convective flux (also referred to as the enthalpy flux) is by 
definition the quantity: 

F c = u (pC p AT) (78) 

where c p = (ds/dh\T) p . Finally, from the definition of Eq. ([71)1 and the set of 
Eqs. (|76 p -([78 p . one derives Fkin oc g A/T F c and, since A oc T/g, we show finally 
that Fkin oc F c . 

In the region where the driving is the most efficient, the total energy flux, 
F tot , is no longer totally transported by convection (that is F c (F tot ). However, 
in order to derive an expression that depends only on the surface parameters of 
the star, we will assume that all of the energy is transported by convection ; that 
is F c « Ftot = crT G 4 ff oc g (L/M) where a is the Steffan's constant. Accordingly, 
Eq. (f75)) can be further simplified as: 

V (x u; 2 8C T c % T 3 g^ « J osc T 7 cS g^ , (79) 



5 for the sake of simplicity we assume here an isotropic medium, accordingly the flux 
of kinetic energy is the same in any direction 
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where we have assumed T = T c g. 

Let now defines v max — lu^ x /2tt the peak frequency associated with V . This 
characteristic frequency can be estimated according to: 

fmax ~ U /A (80) 

where the quantity uo/A is estimated in the layer where uq is maximum. Using 
similar simplifications as used previously for V, we can show that 

^max«3(rclf/p) 1/3 - (81) 

where p is the mean density at the photosphere. We assume that p is equal to 
the star mean density, that is p w M/R 3 oc g/R. Accordingly, we then derive 
from Eq. (EJJ) and Eq. (JHTJ : 

V^oc^) 23 ' 12 9 - 3 MV 3 , (82) 

where M is the stellar mass. For main sequence stars lying in the domain where 
solar-like oscillations are expected, M 1 / 3 varies very slowly such that it can be 
ignored in Eq. (|82|) . Then, Eq. (|82|) can finally be simplified as: 

P raax oc (T c 4 ff ) 2 <T 3 . (83) 

We now clearly see from Eq. (|83|) that "P m ax as expected increases with increasing 
Ftot — o T c 4 ff and decreases with increasing g. 



10.2 Mode surface velocity 



Prior to the CoRoT mission, only crude and indirect derivatio ns of the aver 



aged mode line width had been proposed for a few stars (see iKieldsen et al 



20051: iFletcher et ail 120061: IKieldsen et""all l2008fh However, for the majority of 
solar-like pulsators observed so far from the ground in Doppler velocity, such 
measurements are not available, only the maximum of the mode surface velocity 
(U max hereafter) is in general accessible. For the numerous solar-like pulsators 
observed from the ground, we must compute the mode surface velocity according 
to Eq. (fT3"|) , which requires the knowledge of not only V but also of the mode 
da mping rates (77 = 7r F) . 

Houdek etafl Jl 999) have computed rj for a large set of main sequence mod- 



els. Using Balmforthl ( 1992a Vs formulation of stochastic excitation, they have 



also computed the mode excitation rates (V). From their theoretical compu- 
tations of V and r — rj/ir, they have derived v s according to Eq. (fT3")) . Their 
theoretical calculations for U max result in a scaling law of the form (L/M)P with 
a exponent (3=1.5 (see Irloudek et al. I ll999h . 

We have plotted in Fig. [Tl] the quantity U max associated with the solar-like 
pulsators observed so far in Doppler velocity. Clearly, U ma x increases as (L/M) 13 
where the expone nt ~ 0.7. A similar s caling law with the exponent 0=1 was 
earlier derived by IKieldsen fc Beddind (l995) from the theoretical calculations 
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bv lChristensen-Dalsgaard k. Frandsen (|l983h . iHoudek et al. I <|l999h 's scaling law 
significantly over-estimates the mode amplitudes in F-type stars. For instance 
for Procyon (T cff ~ 6480 K, L ~ 6.9 L© and L/M ~ 4.6), this scaling law 
over-estimates V max by a factor ~ 4. 




Lorentzian ^ k 
Goussion x k 



3 = 0.7) 
3=1.0) 



(L/M) [solor units] 



Fig. 11. Ratio between V ma x the maximum of the mo de velocity relative to the observed 
solar value (V^ ) ax = 25.2 cm/s for 1=1 modes, see lKieldsen et all l2008h . Filled dots 
correspond to the stars for which solar-like o scillations have been detect ed in Doppler 
velocity (see a detailed list of references in iBedding fc Kieldsenl 120071 ) . The lines - 
except the dot-dashed line - correspond to the power laws obtained from the predicted 
scaling laws for "Pmax and estimated values of the damping rates 7/ ma x (see text for 
details). Results for two different eddy time-correlation functions, \k, are presented: 
Lorentzian Xfc (solid line) and Gaussian \k (dashed line). 



Samadi et alJ (12007 ) have derived Vm ax using mode damping rates computed 
bv lHoudek et al.l(|l999 ) and the different scaling laws found for 'Pmax oc (L/M) a . 
They also found that Vmax scales as (L/M) 13 . This is not surprising since 'Pmax 
varies as (L/M) a . Furthermore, the exponent (3 is found to depend significantly 
on the choice of Xk- P = 0.7 for a Lorentzian \k and (3 = 1 for a Gaussian \k- 
As shown in Fig[TlJ the best agreement with the observations is found when a 
Lorentzian \k is assumed. On the other hand, assuming a Gaussian Xk results 
in a larger exponent (3. When theoretical mode amplitudes are calibrated with 
respect to the solar mode amplitudes, calculations based on a Gaussian \k over- 
estimate the amplitudes of solar-like pulsators significantly more luminous than 
the Sun. 

Theoretic al calculations by Hqudek et al. ( 19991 ) assume a Gaussian \k- Then 
according to ISamadi et al.1 (|2007l )'s results, the too large value found for (3 by 



34 R. Samadi 



Houdek et al. ( 1999h c an partially be explained by the use of a Gaussian Xk- 



However, according to iHoudek their too high value of (5 might be ex- 

plained essentially by the mode damping rates that could be under-estimated by 
a factor ~ 1.8. 



11 Discussion and perspectives 



The way mode excitation by turbulent convection is modeled is still very sim- 
plified. As discussed below, several approximations must be improved, some 
assumptions or hypothesis must be removed. 

As seen in Sect. 15. 2\ the driving efficiency crucially depends on the eddy 
time-correlation (xfe)- Current models assume that Xk varies with u> in the same 
way at any length scale. At the length scale of the energy bearing eddy, there are 
some strong indications that Xk is Lorentzian rather than Gaussian. However, 
at smaller scale, it is not yet clear what is the correct description for Xk- Use 
of more realistic 3D simulations would be very helpful to represent the correct 
dynamic behavior of the small-scales. 

Current theoretical models that include the entropy fluctuations in the driv- 
ing assume that the entropy fluctuations behave as a passive scalar (see Sect. l4.3[) . 
As a consequence, cross terms between Sr and Ss vanish. This is a strong hy- 
pothesis that is unlikely to be valid in the super-adiabatic part of the convective 
zone where driving by the entropy is important. Indeed, the super-adiabatic 
layer is a place where the radiative losses of the eddies are important because of 
the optically thin layers. Assuming that the entropy (or equivalently the tem- 
perature) is diffusive (Eq. (J33J ) is no longer valid. Furthermore, departure from 
incompressible turbulence is the largest in that layer and, accordingly, the cross 
terms between Sr and Ss no longer vanish (see SG). Therefore, the passive scalar 
assumption is not valid in the super-adiabatic layers. To avoid this assumption, 
one needs to include the radiative losses in the modeling. 

One other approximation concerns the spatial separation between the modes 
and the contributing eddies. This approximation is less valid in the super- 
adiabatic region where the turbulent Mach number is no longer small, in par- 
ticular for high i order modes. This spatial separation can however be avoided 
if the kinetic energy spectrum associated with the turbulent elements (E(k)) is 
properly coupled with the spatial dependence of the modes (work in progress). 



The CoRoT mission, launched 27 December, is precise en ough to detect solar - 
like oscillations with amplitudes as low as the solar p modes ([Michel et al.ll2008l ). 
Furthermore, thanks to its long term (up to 150 days) and continuous obser- 
vations, it is possible with CoRoT to resolve solar-like oscillations, and hence 
to measu re not only the mode amp litudes but also directly the mode linewidths 
(see e.g. Appourchaux et al. 2008). Similarly as in the case of the Sun, it is 
now possible with CoRoT to derive direct constraints on V for stars with differ- 
ent characteristics: evolutionary status, effective temperature, gravity, chemical 
composition, magnetic field, rotation, surface convection, ... etc. We emphasize 
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below some physical processes and conditions that we expect to address thanks 
to the CoRoT data. 

Some solar-like pulsators are young stars that show rather strong activity 
(e.g. HD 49933, HD 181420,HD 175726, HD 181906, ...). A high level of ac- 
tivity is often linked to the presence of strong magnetic field. Effects of the 
magnetic field are not taken into account in the calculation of the mode excita- 
tion rate s. A strong magnetic fi e ld can more or less inhibit convective transport 
see e.g. Proctor fc Weisslll982t IVogler et al.ll2005l ). Furthermore, as shown by 



Jacoutot et al. ( 2008a[ ). a strong magnetic field can significantly change the way 



turbulent kinetic energy is spatially distributed and leads to a less efficient driv- 
ing of the acoustic modes. In that framework, the CoRoT target HD 175726 is 
probably an interesting case. Indeed, this star shows both a particularly high 
level of activit y and solar-like osc illations with amplitudes significantly lower 
than expected ( Mosser et al. 20091 ). 

Youn g and active stars rota te usually faster than the Sun. As shown re- 
cently bv lBelkacem et al.l (2009aJ), the presence of rotation introduces additional 
sources of driving. However, in the case of a moderate rotator such as HD 49933, 
these additional sources of driving remain negligible compared to the Reynolds 
stress and the entropy source term. On the other hand, the presence of rotation 
has an indirect effect on mode driving through the modification of the mode 
eigenfunctions. An open issue is: will the CoRoT or the Kepler m ission be able 
to test the expected effect of rotation (see iBelkacem et al]|2009al) ? 

Solar- like oscillations have now been firmly detected in se veral red giant stars, 
from both Doppler velocity measurements (see the review by iBedding fc Kicldscn 



2006) as well as from spa ce based photometry measurements (jBarban et alJ 



2007; de Ridder et al.ll2006l ). More recently, detection of solar-lik e oscillations by 



CoRoT in a huge number of red giant stars has been announced bv lde Ridder et al 



(2009). Why look at solar-like oscillations in red giant stars? Toward the end of 
their lives, stars like the Sun greatly expand to become giant stars. A consequence 
of this great expand, is the existence of a very dilute convective envelope. A low 
density favors a vigorous convection, hence higher Mach numbers (M t ). The 
theoretical models of stochastic excitation are strictly valid in a medium where 
Mt is - as in the Sun and a Cen A - rather small. Hence, the higher M t , the 
more questionable the different approximations and the assumptions involved in 
the theory. Hence, red giant stars allow us to test the theory of mode driving by 
turbulent in more extreme conditions. 

Finally, most of theories of sto c hastic excitation a re developed for radia l 
modes only. iDolginov fc Muslimovl (|l984f ). GMK and IBelkacem et all (|2008l ) 
have considered the non-radial case. There are interesting appl ications of such 



non-ra dial formalisms, for instance the case of solar g modes (IBelkacem et al 
2009bh . but also g modes in massive stars that ca n in principle be excited in 



their central convective zones (jSamadi et al.ll2009ah . 
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